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ABSTRACT 

In previous warm-up chemical models of the low-mass star-forming region 
L1527, we investigated the evolution of carbon-chain unsaturated hydrocarbon 
species when the envelope temperature is slightly elevated to T w 30 K. These 
models demonstrated that enhanced abundances of such species can be explained 
by gas-phase ion-molecule chemistry following the partial sublimation of methane 
from grain surfaces. We also concluded that the abundances of hydrocarbon 
radicals such as the C n H family should be further enhanced as the temperatures 
increase to higher values, but this conclusion stood in contrast with the lack of un- 
ambiguous detection of these species toward hot core and corino sources. Mean- 
while, observational surveys have identified C2H, C4H, CH3CCH, and CH3OH 
toward hot corinos (especially IRAS 16293-2422) as well as towards L1527, with 
lower abundances for the carbon chain radicals and higher abundances for the 
other two species toward the hot corinos. In addition, the Herschel Space Tele- 
scope has detected the bare linear chain C3 in 50 K material surrounding young 
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high-mass stellar objects. To understand these new results, we revisit previ- 
ous warm-up models with an augmented gas-grain network that incorporated 
reactions from a gas-phase network that was constructed for use with increased 
temperature up to 800 K. Some of the newly adopted reactions between carbon- 
chain species and abundant H2 possess chemical activation energy barriers. The 
revised model results now better reproduce the observed abundances of unsatu- 
rated carbon chains under hot-corino (100 K) conditions and make predictions for 
the abundances of bare carbon chains in the 50 K regions observed by Herschel 
HIFI. 

Subject headings: Astrochemistry, ISM: molecules, ISM:individual (L1527, B228, 
IRAS 16293-2422, W31C, W49N), stars:formation 



1. Introduction 



The synthesis and evolution of carbon-chain species have once again become topics 
of interest, specifically pertaining to early stages of low-mass star formation in which the 
envelope temperature is warmer than the parent dark cloud temperature of 10 K. Recent 
observations indicate that two specific envelope stages for low-mass protostars, hot and 
lukewarm corinos, show differences in their abundances of carbon chain species. 

Hot corinos are relatively condensed (tih ~ 10 7 cm -3 ) and hot (T > 100 K) re- 
gions, and feature a molecular composition rich in hydrogen-rich organic species as well 
as oxygen-bearing species. They can be considered low-mass analogs to the hot cores that 
surround high-mass star-forming regions. Some of t he first hot corinos to be categorized 



were IRAS 16293-2422, NGC 1333 IRAS 4A and 4B flBottinelli et al.l 12004 IjQrgensen et al. 



20041 ; iMaret et al.l 120041 ; ICazaux et al.l 120031 ) . As in hot core s, the complex composi tion is 
the combined result of grain surface and gas phase processes (IGarrod fc Herbstll2006l ). 



Observations of L1527, an envelope surrounding the transitional Class 0/1 protostar 
IRAS 04368+2557 in Taurus, indicated that enhanced abundances of unsaturated (hydrogen- 
poor) hydrocarbons are likely related to a slightly elevated envelope t emperature of T ^ 
30 K via a chemistry labeled "warm carbon-chain chemistry" (WCCC) (jSakai et al.ll2008bl ). 
Although the chemical composition resembles that of a dark cloud core, the composition of 
WCCC regions can also be attributed to a process beginning with grain surface chemistry. 
The elevated temperature allows methane to evaporate from the ice mantles surrounding 
dust particles and act as a precursor for a carbon-chain rich ion-molecule chemistry. The 
resulting chemical composition is a defining characteristic of the source, which we designated 
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a "lukewarm corino" . The idea of WCCC has proved somewhat controversial, however, and 
alternative explanations, such as standard gas-phase ion- molecule chemistry in a collapsing 
source, have been advanced (jCharnley fc Cordinerl |2010| ). A recent interferometric study 
provides evidence for WCCC chemistry within the central protostar and normal ion-molecule 
chemistry in cold cloud remnants 



(ISakai et al.l l2010h . It has not been firmly established if the lukewarm corino phase 
is one step on an evolutionary track leading to the hot corino phase, or if it represents a 
terminal stage. 

Since the observation of L1527, additional observations have identified a similar kinetic 
temperature and similar enhancements for hydrocarb on chains towards the source IRAS 
15398-3359 in Lupus, alternatively identified as B228 (|Sakai et al.ll2009af) . so that B228 is 
the second lukewarm corino to be studied. Also, ISakai et al.l ( I2009al ) reported detections of 
some carbon chains, especially C4H, toward a number of hot corinos at temperatures near 100 
K, but with abundances significantly lower than those in the two lukewarm corino sources, 
suggesting a diminution in abundance with increasing temperature. More recently, linear C3 
has been detected in absorption in regions at a te mperature of ~ 50 K a gainst continua in 
the THz region from the sources W31C and W49N (IMookerjea et al.ll2010l . hereafter referred 
to as the HIFI sources). 

We previously reported the use of a gas-grain chemistry network with a simple warm-up 
following a cold stage to si mulate the evolutio n of molecular abundances in L1527 as a result 
of increasing temperature (IHassel et al.ll2008l)f hereafter H HG) . A comparison of m odel abun- 
dances with the observations of Ijprgensen et al.l (120041 ) and ISakai et al.l ( 12008al Jbl) showed 
that the best overall agreement occurs at times when the temperature approaches T = 30 
K, the indicated temperature for the envelope. Reasonable agreement with observation in 
L1527 was found near this temperature whether or not the warm-up was terminated at 30 K 
or continued to higher temperatures. Similar chemistry was al so found in the hydro dynamic- 
chemical model of prestellar to protostellar core formation of lAikawa et al.l (120081 ) . 



Some of the warm-up models were extended to maximum temperatures of 100 and 
200 K, in a similar fashion to earlier hot core models, in which it was shown that the 
standard oxygen-containing complex mol ecules in such sources could be partially prod uced by 
surface chemistry followed by desorption (IGarrod fc Herbstll2006t iGarrod et al.ll2008f ) . These 
prior studies emphasized the typical hot-core species instead of the calculated abundances 
of carbon-chain species. It was later noted that in all of our warm-up models, carbon chain 
species are quite abundant at the higher temperatures for long periods of time. This result 
contrasts with the lack of strong detections of such species toward hot co re and hot corino 
sources, with the general exception of HC3N (IHerbst fc van Dishoeckll2009[ ). In the survey of 
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Sakai et al.l (j2009al ). the fractional abundance values obtained for C4H vary, but are generally 



smaller than X Ci}1 = Nc 4 u/^h = 10" in hot corinos, as compared with X C4H ~ 2-3 xlO 
for the two lukewarm corino sources. The observed hot corino result is smaller than the peak 
fractional abundance of Xc 4 h > 10~ 8 found in the extended warm-up models of HHG. 

Some other observational studies have been directed at the question of unsaturated 
molecules in star-forming reg ions. In their review of the chemistry of complex molecules, 
Herbst fc van Dishoeckl J2009I ) mentioned a study of the C4H/CH3OH and HC5N/CH3OH 
abundance ratios in seven hot cores that resulted in upper limits of 0.01-0.1 for the ratios (S. 
and J. J0rgensen-Bisschop, private communication). In addit ion, an average abun dance of 
-/Vch 3 oh/^Vh = 9.5 x 1CT 7 was reported for the seven hot cores (IBisschop et al.ll2007l ). leading 
to upper limits for C 4 H and HC 5 N high enough to overlap with the T max = 100 & 200 K 
peak results of HHG for both species. 

It would appear at present that C4H is detected in largest abundance in the two known 
lukewarm corinos, that it is detectable in hot corino sources although to much lower abun- 
dance than predicte d, and that it has a high upper limit with respect to methanol in hot 
cores. The survey of ISakai et al.l (l2009bl ) also included the detection of CH 3 OH toward LI 527 
with a fractional abundance of ATch 3 oh ~ 10" 9 , similar to what is observed in cold cores 
such as TMC-1. Unlike C4H, and possibly other radicals of this class, methanol is a molecule 
typically associated with high gas-phase abundances only in hot cores and corinos, where 
it is desorbed from grains. Although the picture is clearly incomplete at the present time, 
comparison of carbon chain species and species known to prefer high temperatures provides 
an initial basis for chemical distinction of lukewarm and hot corino sources. A lthough one 
can thi nk of a lukewarm corino as a prior stage of evolution than a hot corino, ISakai et al. 
( 12009bl ) have also suggested that the difference between lukewarm and hot corinos has to 
do at least partially with the rate of collapse of the protostar; the slower the rate, the more 
chance for CO to form in the gas, accrete onto grains, and lead to oxygen-containing species; 
the faster the rate, the greater the chance that C atoms land on grains and form methane. 



In a recent study, lHarada et al.l ( 120101 ) reported the extension of the temperature range 
of the OSU gas-phase chemical model to 800 K by the inclusion of many new reactions, 
including those with non-zero activation barriers. Select reactions might also be relevant 
at hot core and corino temperatures. In particular, reactions between the carbon chain 
radicals C n H and C n and H 2 were omitted from previous networks because they possess 
barriers of Ex > 10 3 K, and are inefficient in cold models. The effect of adding these gas- 
phase reactions to our warm-up gas-grain model network would likely be to decrease the 
calculated abundances of these radicals at higher temperatures, so as to make the results 
more in line with hot corino observations of C4H. On the other hand, activation energy 
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barriers for reactions between normal molecules in ground electronic singlet states, such as 
cyanopolyynes, and H 2 should be much larger so that these reactions are still not likely to be 
important in hot corinos. The updated reaction network should not affect the abundances 
at lukewarm corino conditions. In this paper, we present revised gas-grain warm-up models 
inc luding addition of g as-phase radical-H2 reactions as well as many other reactions added 
by lHarada et al.l (120101 ). and use the results to determine a revised temperature dependence 
of the abundances of carbon-chain species in the range from 30 K (lukewarm corinos) to 
100 K (hot corinos). 



2. Model 



We utilize an updated version of the Ohio State University (OSU) gas-grain code to 
follow the evolution of abundances duri ng a warm-up to simu l ate the central regions sur- 
rounding a burgeoning protostellar core (IHasegawa et al.lll992l ; iGarrod fc Herbstll2006l ). In 
the model, a one-point parcel of material 

undergoes a gradual temperature increase following an initially cold period of To = 10 K 
with a duration of t = 10 5 yr. This initially cold period allows for the accumulation of 
significant solid-phase abundances of H 2 0(s), CO(s) and CH 4 (s). Varying the length of 
the cold period between 10 4 -10 6 yr did not significantly change the results of HHG, so we 
presently consider models that begin warming up at 10 5 yr in the interest of making direct 
comparisons with prior results. 

We consider constant gas densities of nn = 2 x 10 5 and 1 x 10 6 cm -3 based on estimates 
for the HIFI and lukewarm corino sourc es, respectively. The same physical parameters and 
initial chemical abundances as found in IGarrod et al.l ( 120071 ) and Tables 1 & 2 of HHG are 
adopted. Our approach uses homogeneous physical conditions at any given time, which i s 
a simplification not found in more complex hydrodynamic approaches (lAikawa et al.l 120081 ). 
The advantage of this approach is that it allows a more detailed look at the chemical processes 
and the roles of individual reactions. 



The warm-up rate is again treated according to 

T = Tq + (T max — Tq) 



T 



where we use values for the maximum temperature T max of 30 K, 50 K, 100 K, and 200 K. 
The 30 K value is used as a terminal temperature for lukewarm corino models, while the 100 
K and 200 K values are used to represent sequences in which a lukewarm corino proceeds 
towards a hot corino. The maximum temperature of 50 K is used for HIFI models. 
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A heating timescal e of th = 2 x 10 5 yr is adopted, with an order of n = 2 following 
Garrod fc Herbstl (120061 ). In addition to the warm-up models, we include a 10 K constant 
temperature model. 

The updated gas-grain network contains 7166 reactions involving a total of 668 gaseous 
and surface species. Among the new gas-phase reactions, the most important are several 
sets of hydrocarbon radical-H 2 gas-phase reactions with moderate barriers, and reactions 
that distinguish the linear/carbene isomers HC3 and H2C3 from the cyclic forms C-C3H and 
C-C3H2. The three new reaction classes that are most important for simulating radical chains 
in hot corinos are listed in Table [TJ The table also includes parameters for computing the 
modified Arrhenius rate coefficient of a gas phase reaction: 



a 



T 



300 K 



cxp 



-Ea 
T 



(2) 



where Ex (more commonly referred to as 7) represents the activation barrier. The rate 
coefficients used for th e first class (CrM + H?) are based o n assorted studies of the exothermic 
reaction with n = 2 (IFarhat et al.l Il993t iHerbstl Il994l for example). Since the reported 
values for the rate coefficient show some dispersion, we have chosen to adopt the modified 
Arrhenius parameters from shock models (IHarada et al.ll2010l . G. Pineau des Forets, private 
communication), which should be reasonably accurate in the temperature range of interest. 
Similarly, the rate coefficients used for t he second class (C n H + H?) a re based on the value 
for n=2 from shock- model compilations fjPineau des Forets et al.lll988l ). This reaction is also 
exothermic. For the third set o f reactions (C n N + H ?), we used the measured values for the 
exothermic CN + H2 reaction (IWoodall et al.l 120071 ). The backwards endothermic reactions 
for these three reaction classes are ignored. 



3. Results 

In this section, we report the effects of the expanded reaction network on various warm- 
up models for different classes of carbon-chain species. We contrast the evolution of these 
species with that of methanol, a species more typically associated with hot corinos and 
hot cores. We consider a broader comparison with observational results for corinos, both 
lukewarm and hot, in §H 
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3.1. Polyyne Radicals: C n H 

Figure [1] shows the evolution of the polyyne radical C4H in assorted models with the 
updated reaction network (solid lines) and the previous results (dashed lines). Panel a) 
contains our constant-temperature, cold (10 K) model, whereas the three other panels contain 
the results of warm-up models, in which the warm-up phase occurs from 1 — 3 x 10 5 yr, 
followed by a phase at the asymptotic temperature of 30 K, 100 K, and 200 K, in panels 
b), c), and d) respectively. The observational abundances for the two lukewarm corinos are 
shown as horizontal lines. The results in panel a) show comparatively low abundances of 
C4H compared with observation for both old and revised models. Panel b), in which the 
asymptotic temperature is 30 K, can be construed as a model of a lukewarm corino that 
remains at 30 K rather than continuing to higher temperatures. Here one can see little 
change between the old and new results, indicating that up to 30 K, the new reactions have 
little effect. Panels c) and d) show our warm-up models leading to T max = 100 K and 200 
K, roughly corresponding to the conditions of hot corinos and hot cores, respectively. The 
C4H abundances rise to values comparable with the peak of the T max = 30 K model as the 
temperatures approach 30 K, in a similar manner to our prior model. As the temperature 
continues to rise, however, a difference between the old and new models emerges. In the 
older calculations, the C4H abundance continues to rise to a peak value that exceeds 10 -8 
and remains larger than 10~ 9 for a time period of t ~ 10 6 yr, whereas in the updated ones, 
the C4H abundance almost immediately declines. 

Figure |2] amplifies the abundance change of C4H and other species (C2H, CH3CCH, 
CH3OH) for the model in which the warm up continues all the way to 200 K. The results for 
C4H appear in panel (b), where some more details of the differences between the previous 
and newer models can be seen in the closer view. After the temperature reaches ~ 80 K, the 
major distinction becomes apparent, as the reaction C4H + H2 — >■ C4H2 + H becomes the 
dominant destruction mechanism in the newer model, and the abundance of C4H declines 
precipitously. 

The abundances for the C n H chains evolve in a similar manner to C4H in the assorted 
models. The constant T = 10 K model results agree between newer and previous models 
to the same degree as the results appearing in Figure and the T max = 30 K warm-up 
model results rise to a similar peak value for a similar duration with bold old and new 
models. Likewise, the T max = 100 and 200 K model results rise to similar peak values as 
the corresponding 30 K results, but differences appear beyond this time in these models as 
a consequence of the inclusion of radical-H2 reactions, which reduce the radical abundances 
quickly in the new calculations. For the case of C3H, there is a distinction in the new 
calculations between the linear (HC3) and cyclic (C-C3H) forms, with the latter having a 
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somewhat higher abundance, as shown in Figure [31 which contains plots of the abundances 
vs time for three warm-up models. As expected, we note a relatively decreased abundance 
at higher temperatures for both isomeric forms as compared with the previous results, as a 
result of the onset of high-temperature reaction with H 2 . 

3.2. Polyacetylenes, Rings, and Carbenes: HC n H, c-C n H 2 , and H 2 C n 

The products of the reactions in Table [1] imply that if these destruction reactions re- 
duce the abundance of C n Ii chains at higher temperatures, there should be a corresponding 
increase in C n H 2 species. This effect is seen at certain times for selected molecules, but 
the C n H 2 species can exist as different isomers, including the non-polar polyacetylene form 
HC n H, the carbene form H 2 C n , and various ring forms c-C n H 2 . This can make verification 
of models difficult both because reaction networks do not always distinguish the isomers 
and because only some of them are polar and have detectable rotational spectra. In our 
previous analysis, we assumed that the carbene forms represent 2% of the abundance of the 
undifferentiated model abundances, and we generally found agreement within an order of 
magnitude between observations and model results with this modification. The assumption 
was based on the assumed rapid destruction with atomic oxygen for the carbene but not the 
more stable acetylenic species. 

We show the calculated evolution in Figure H] for c-CsH 2 and H 2 C3 for various warm- 
up models in the newer network and for the previously undistinguished species CsH 2 from 
HHG. The results are more complex than expected. The more stable c-C3H 2 isomer is more 
abundant than the carbene isomer at most times in the assorted warm-up models, although 
often by less than the factor of 50 assumed in HHG. Unexpectedly, the cyclic isomer can be 
more abundant than the undifferentiated species in the HHG model, which is found to be 
a secondary effect caused by a number of new barrierless reactions in the updated network. 
In addition, at late times in the T max = 30 K model an inversion of the H 2 C3 and c-C3H 2 
abundances is apparent, an effect also caused by new reactions other than those with barriers. 

The initial evolution during the warm-up progresses in a similar manner in panels b) 
and c) for the T max = 100 K and 200 K models, but an increase in the peak abundances of 
the species becomes apparent as a result of desorption from the grain surfaces. Differences 
in surface rates prior to the onset of desoprtion result in more C3H 2 (s) than H 2 C3(s), and 
ultimately affect the gas phase abundances during desorption. The difference in the abun- 
dance of CsH 2 between current and prior results can be attributed to the onset of C3H + H 2 
— > CsH 2 + H for T > 80 K. Panel c) illustrates a further difference between CsH 2 and H 2 C3 
in the current model as the destruction rate H 2 C3 + H 2 — > C3H3 + H becomes a dominant 
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process as T > 140 K. The C3H3 then reacts to preferentially form C3H2 rather than H2C3. 
The energy barrier for this rate of = 1500 K is significantly smaller than the value of 
27,100 K for the corresponding destruction reaction involving C3H2. This difference has the 
overall effect of converting not only C3H but also H2C3 into C3H 2 at high temperatures. 

The species C4H2 and C6H2 also exist in different isomeric forms, but these were not 
differentiated in the updated reaction set. These species differ from C 3 H 2 because only the 
carbene form has been observed. A comparison of model results shows little discernible 
difference between the current results and previous HHG models for C4H2 at lower tem- 
peratures, but some difference at late times in the T max = 200 K model. The new CeH 2 
abundances show a more distinct enhancement over the previous ones during warm-up as 
T > 80 K. This suggests that the updated reaction set drives the longer-chain C n H material 
toward C n H 2 species at high temperature conditions. Given the comparison of undifferenti- 
ated model results and observations for the H2C4 and H2C6 isomers, it is likely that much 
of the carbon-chain material is deposited into the unobservable forms HC n H. 



3.3. Cyanopolyynes: HC n N 

The absence of viable destructive reactions with H 2 for cyanopolyyne species implies 
that these should not be depleted at high temperatures. Instead, the reactions in Table [T] 
suggest that radical C n N chains should be converted to HC n N species as T > 80 K. This 
further suggests that the C n N evolution might be similar to that of C n H, while the HC n N 
evolution might resemble C n H 2 . 

The warm-up model results for C5N and HC5N are illustrated in Figure O and compared 
with previous results from the HHG analysis. In panel a), it can be seen that the old and 
new results for the T max = 30 K models differ for both species, and to a larger degree than 
for C 4 H. The higher temperature results in panel c) do support the conjecture that C 5 N will 
peak, and then rapidly decline with increasing T in favor of HC5N, in stark contrast with 
previous results. This trend appears also in panel b), but to a lesser extent. 



3.4. Methanol: CH 3 OH 



In contrast with the unsaturated carb on chain species, CH3OH is more ty pically asso- 
ciated with hot core/corino environments (ILovas et al.lll982t iBlake et al.l 119871 ). Thus, the 
abundance is anticipated to increase with increasing temperature, and should not be read- 
ily destroyed by the inclusion of the high temperature network. We include a comparison 
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between the old and new results for the T max = 200 K models in Figure [2H. Some differ- 
ence in abundance can be seen between the models for the range T ?s 20 — 80 K, which 
can be attributed to the inclusion of a photodesorption process in the updated network 
(lOberg et al.ll2009l ). As the temperature rises beyond this range, the abundance of CH 3 OH 
rises by nearly four orders of magnitude as a result of thermal desorption from grain surfaces. 
Little discernible difference is noted between the two models until later times following the 
warmup. 



3.5. Bare Radical Carbon Chains: C n 

Our previous analysis included but did not focus on bare carbon chains, C n , pri- 
marily because detecti ons of such species had not been made toward L1527. Since then, 



Mookerjea et al.l (120101 ) have observed vibrational-rotational transitions of C3 in absorption 
against hot-core type regions W31C and W49N using the HIFI detector on Herschel. The ab- 
sorption is thought to arise from a cooler envelope surrounding the hot cores at an estimated 
temperature of ~ 50 K. Their results followed an earlier, more tentative interstellar detection 
towards Sgr B2 at low resolution. While no larger bare carbon chains have been detected 
by Herschel, we have re-computed results to n = 10 and show results between n = 3 — 6 in 
Figure E]for a m odel at a density »h = 2 x 10 5 cm -3 and a warm-up to T max = 50 K based on 



the estimates of iMookerjea et al.l (120101 ). There is not a significant difference between the C n 
abundances from the new models and our prior results because the T max = 50 K is not hot 
enough to make the C n + H2 — > C n H + H reactions relevant as a destruction mechanism. 
Therefore, we switch focus and consider the role of surface chemistry in their formation. 

For each bare carbon chain, we plot the new results in Figure E] including surface chem- 
istry and results for an analogous model without surface chemistry. We note that with the 
gas-grain models, the pre-warm-up abundances of the odd-numbered chains C3 and C5 ex- 
ceed the secondary peaks caused by the warm-up, which starts at t = 10 5 yr. In the purely 
gas-phase models, on the other hand, there is no warm-up peak at all; rather, there is a 
precipitous drop in the abundances of the bare carbon chains. These results illustrate that 
the warm-up peaks clearly require surface chemistry. However, as a more detailed analysis 
shows, the surface chemistry is not direct, but serves mainly to produce methane during 
the cold era, which evaporates as the temperature nears 30 K, an d acts as a precur sor for 



another round of ion-molecule chemistry leading to carbon chains (jSakai et al.ll2008bf ). Like 
the other carbon-chain radicals, the primary re-generation mechanism for the C3 warm-up 
peak is directly related to acetylene derived from methane, and occurs via the gas-phase 
reaction C+C2H2 — >■ C3 + H2. During the warm-up to 50 K, all four chains reach an abun- 
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dance peak of X ~ 10 8 for approximately the same time frame of several times 10 5 yr, with 
only a gradual decline in the peak height as the chain length increases. 



4. Comparison with Observations 

The motivation for including high-temperature reactions stems partially from an in- 
crease in relevant observational data. This increase includes new detections toward L1527, 
hot corinos, and HIFI sources, as well as a new candidate lukewarm corino. The expanded 
data set allows for the compositional comparison of these sources. 



4.1. Lukewarm Corinos 



Observed fractional abundances with respect to % for the lukewarm corino sources 
L1527 and B228 are listed in Table [2j In addition to the new observations of B228, this table 
includ es some new and revisited observations toward LI 52 7 (IHirota et al.ll2009t ISakai et al. 
2009aU bl). Some updates include the distinction between the c-C 3 H 2 and H 2 C 3 isomers, 
and the detections of I-C3H and CH3OH. To convert observations from column densities to 
fractional abundances we used a total column density of Nu = 6 x 10 22 cm -2 for LI 527 



( jjprgensen et al.l 120021 ). Based on dust continuum flux measurements, ISakai et al.l ( 12009al ) 
adopted the same value for B228. 

The mean confidence level, k, was used to quant ify the agreement be tween model results 
and observed abundances based on the approach of iGarrod et al.l (120071 ). With their defini- 
tion, an average confidence level of unity refers to perfect agreement, an average confidence 
level of 0.317 refers to an average deviation between observed and calculated abundances 
of an order of magnitude, while an average confidence level of 0.046 refers to an average 
deviation of two orders of magnitude. This mean confidence level is calculated as a function 
of time for our models; the best value is referred to as ft max , and it occurs at a time t opt and 
corresponding temperature T. We also determine the number of species that agree within 
an order of magnitude at that time. We used our 10 K constant-temperature model and 
three warm-up models (T max = 30 K, 100 K, and 200 K) to study L1527. The results appear 
in Table H 

The inclusion of new reactions and new detections into the comparison has somewhat 
decreased the overall level of agreement from our previous work. Part of the problem is 
the need to fit the abundances of the carbene structures without the simple approximation 
that they have 2% of the polyacetylene abundances. First, we find that the t opt values, 
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and thus the associated temperatures, do not change for the three warm-up models when 
compared with previous results. The values of ft max for the T max = 30 K, 100 K, and 200 
K models decrease slightly in the new analysis to 0.599, 0.580, and 0.577, corresponding 
to overall agreement within a factor of 3.5-4 at optimum temperatures of 25 K, 26 K, and 
27 K, respectively. In terms of number of molecules, the warm up models fit 19, 20, and 
19 species out of 26 to within an order of magnitude. This level of agreement is not much 
better than the agreement with the 10 K constant-temperature model of K max = 0.514, which 
occurs at a time of 4.2 x 10 4 yr, slightly before warm-up starts. Here 19 species are fit to 
within an order of magnitude. This model, however, still significantly under-produces some 
major characteristic carbon chain species: C2H, C4H, HC5N, HC7N, and HC9N at the t opt 
value, and must therefore be regarded as somewhat inferior for these species, whatever its 
reliability for non-carbon-chain species. 

The newly observed source B228 is presently limited to seven species, each of which 
agrees in abundance with the corresponding LI 52 7 values within a factor of 5 or better. The 
new warm-up models fit five of these, but over-produce H2C3 and H2C4 by more than an 
order of magnitude. The constant temperature model is mostly consistent with the L1527 
comparison, but it slightly under-produces CH3CCH toward B228 while remaining at the 
borderline for agreement with L1527. In summary, the updates of new species and reactions 
have somewhat altered the model results for L1527, but have not significantly changed the 
overall conclusion that the fits of warm-up models to L1527 observations near T ~ 25 K are 
a preferable explanation for the observed composition over a constant te mperature m o del at 



T = 10 K. The recent interferometric observations of L1527 conducted by lSakai et al.l (120101 ) 
show that the distributions of C2H and C4H have extended structures as well as an enhanced 
component, which indicate both a regenerative component, fit by our warm-up models, and 
a starless-core phase, corresponding to remnant cold species. 

Based on the fact that the fits to the abundance data L1527 and B228 do not change 
much when the maximum temperature is raised from 30 K to 200 K, it is unclear whether 
lukewarm corinos represent a terminal evolutionary stage at 25 - 30 K or an intermediate 
step of an unfinished process l eading to ho t corin os, or some other stage. It is clear from the 



recent interferometric study of ISakai et al.l (120101 ) that there are depletions of C2H, C4H, and 



C-C3H2 toward the innermost regions of L1527, providing some initial evidence for a third 
chemically distinct region, possibly a hot corino. 
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4.2. Hot Corinos 



In Table [31 we compile observed abundances toward the hot corinos IRAS 16293-2422, 



NGC 1333 IRAS 4A and 4B flSakai et al.ll2009al . 120061 : iBottinelli et al.ll2004l ; Ijorgensen et al. 



20041 ; iMaret et al.ll2004l ; ICazaux et al.ll2003f ) . It can be seen that the first of the three is chem- 
ically richer than the other two, although the three share several chemical characteristics: 
for example, both IRAS 16293-2422 and NGC 1333 IRAS 4B have measurable abundances 
of C4H, which are however two orders of magnitude below that of the lukewarm corinos 
(ISakai et al.ll2009al ; ICazaux et al.ll2003l ). On the other hand, C2H and CH3CCH have only 
been detected in IRAS 16293-2422; while the former has a much lower abundance than seen 
i n the lukewarm co r inos, the more satur ated CH 3 CCH has a considerably higher abundance 
flSakai et al.l l2009at ICazaux et al.ll2003l ). In addition, all three contain high abundances of 
oxygen-containing organic species, such as methanol, and methyl formate. 

In Table El we also present the results of the T max = 100 and 200 K models at T = 50 , 
100, 150 and 200 K to examine the effect of temperature on the evolution of the three 
carbon-chain species and methanol, where we use the results of both models to report a 
range of abundance at 50 and 100 K, and report only the results from the 200 K model 
for T = 150 and 200 K. These results can be compared with those for warm-up models of 
lukewarm corinos, shown in Table [21 where the best results all hover in the temperature 
range 25 — 27 K. Unlike the case of lukewarm corinos, we are not interested in determining 
the best models overall for hot corinos since our gas-grain network does not contain enough 
chemistry to simulate the abundances of the complex oxygen-con taining organi c mole cules 
found in these so urces. This chemistry is studied in the papers of iGarrod et al.l ( 120081 ) and 
Laas et al.l J201l[ ). 



The lukewarm/hot corino distinction is also explored in the solid lines (new model 
results) in Figure El which includes the T max = 200 K models and observed abundances 
for methanol and the three carbon-chain species C2H, C4H, and CH 3 CCH. In comparing 
lukewarm and hot corinos with this warm-up model, we are tacitly assuming that lukewarm 
corinos are earlier and hot corinos later evolutionary stages towards low-mass star formation. 
As the plot illustrates, the radicals C2H and C4H are observed to have higher abundances 
in the lukewarm corinos, while the more hydrogenated CH3CCH and CH3OH have lower 
abundances. The revised model results reproduce this observation. Specifically, the best fit 
to the abundances of the four species in L1527 occurs for T = 30 — 40 K, which is slightly 
higher than the best fit to all species in L1527, which occurs at 25 K as listed in Table El 
On the other hand, the best fit to the four species in the hot corino IRAS 16293-2422 occurs 
at 150 K with some uncertainty, due to sharp temperature variations around this value for 
individual species. We have no overall best fit for this source, but we can conclude from the 
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tabulated results in Table [3] that at T = 100 K all four calculated abundances lie within one 
order of magnitude of their observed values and at higher temperatures and later times, three 
of the four are fit to this criterion. The T = 50 K models still over-produce the polyynes 
and under-produce the hydrogenated species. 



4.3. HIFI Sources 

The HI FI detection of C3 tow ard dense 50 K regions in front of massive star form- 



ing regions (IMookerjea et all |2010| ) may sample gas similar to that in lukewarm corinos. 



In Figure El we includ e the observed abundance of X (C3) ~ 10 8 for W31C and W49N 



( IMookerjea et all 120101 ) along with new warm-up model results with a maximum tempera- 
ture T max = 50 K for C3 - C6 both with and without surface chemistry. While the model 
result without surface reactions crosses the observed abundance for C3 on a steep downward 
trend with increasing time as 50 K is just reached, the model including the surface chem- 
istry agrees with the observation within a factor of 2 or better for a period of approximately 
6 x 10 5 yr, clearly in better agreement with observation. The abundance vs time plots for 
the larger bare carbon chains show an interesting trend in which the species reach a similar 
peak abundance of X (C n ) ~ 10~ 8 within the same approximate time period, with a slight 
decline in peak abundance as chain length increases. This result suggests that bare chains 
up to lengths of 8 or 9 carbons should be abundant enough for detection, provided that their 
spectral vibrational-rotational features fall within the detectable range of an instrument. 
The abundance of C3 is also included in Table [2] at the various t opt values. Although it has 
not yet been reported toward L1527 or B228 and cannot be formally considered to agree, the 
values corresponding to the warm-up models are consistent with Xq 3 = 10~ 8 in these models. 
Additional possible tests to verify this result include searching for bare carbon chains toward 
lukewarm corino sources, or searching for other species in the T ~ 50 K regions surrounding 
the HIFI sources. 



Discussion 



We have calculated gas-grain ch emical models using an expanded reaction network that 
includes non-zero barrier reactions (IHarada et al.l 120101 ) in order to better understand the 
abundances of carbon chains and other species in star-forming regions from lukewarm to hot 
corinos and in dense, cool regions detected by Herschel surrounding young high-mass stellar 
objects. Our models include a cold core period at 10 K, followed by a warm-up period to a 
maximum asymptotic temperature ranging from 30 K to 200 K. The specific temperature of 
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a source can either be associated with an asymptotic temperature or a temperature in the 
range from 10 K to the asymptotic temperature of a model. It is not presently known if a 
lukewarm corino such as L1527 (T « 30 K) represents a terminal stage of chemical evolution 
or a transient phase on an evolutionary track toward the hot corino stage (100-200 K), but 
the present models allow either scenario, because, 

generally speaking, the optimum results we get are not strongly dependent on whether 
a source temperature is considered to be an intermediate temperature or an asymptotic 
temperature. Recent observations show that a range of p hysical conditions and chemical 
compositions can exist not only toward comp lex hot cores (IMookerjea et al.ll2010l ) but also 
towards lukewarm corinos (ISakai et al.l l2010h , and the results provide a basis for future 
modeling investigations. The present models are better suited to address the more specific 
question of whether large abundances of carbon chain species at T > 100 K exist. 



Our four main conclusions follows: 



1. The inclusion of new reactions into our gas-grain network does not significantly affect 
the results and conclusions obtained from previous warm-up model fits to lukewarm 
corino observations. As a benchmark calculation, we revisited the simulations of the 



lukewarm corino L1527 and evaluated fits using the same criteria as lHassel et al.l ( 120081 ). 
The times of peak agreement, t opt , and associated temperatures (25 - 30 K) do not 
change as a result of using an expanded reaction network with reactions containing 
barriers. The values of the optimum mean confidence level, K max , are only slightly 
diminished by comparison with the previously reported values, partly because of the 
addition of more species to the observation set. Thus, the idea that the chemistry 
of a lukewarm corino involves a warm regenera tion of carbon chains resulting from 



CH4 desorption at 25 - 30 K ( ISakai et al.ll2008bl ). incorporated into both our old and 
updated models, is supported. 

2. The inclusion of the new reactions tends to lower carbon-chain radical abundances 
at higher temperatures. Prior models showed significantly increased abundances over 
significant periods of time for several carbon-chain species including the C n H and HC n N 
families above the L1527 observed values as the temperature increased from lukewarm 
corino (T ~ 30 K) toward hot corino/hot core levels (T > 100 K). In the updated 
networks, reactions with activation energy barriers between molecular hydrogen and 
polyyne radicals become the dominant destruction mechanism for these species in star- 
forming regions as temperatures exceed T > 80 K. The inclusion of these reactions 
reduces the abundances of most carbon chain radical species at higher temperatures. 



3. Most importantly, the new observations and models begin to chemically distinguish 
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lukewarm and hot corinos. The species C2H, C4H, CH3CCH, and CH3OH have been 
observed with abundance differences of approximately two orders of magnitude between 
these types of sources, and previous models could not account for the depletion of 
C 2 H and C 4 H at larger temperatures. Updated results at increased temperatures of 
T > 80 K show agreement with the diminished abundances of C2H and C4H observed 
toward hot corinos compared with lukewarm ones. The reactions most responsible 
for diminishing the abundances of the carbon chain (polyyne) radicals C2H and C4H 
and distinguishing lukewarm from hot corinos in their carbon chain inventories are 
those with molecular hydrogen, listed in Table [TJ Even though they possess small 
barriers, these reactions are still competitive at temperatures above 80 K. The model 
results indicate that the "excess" polyynes are primarily converted to unobservable 
polyacetylene forms of C n H 2 , although it is not currently possible to directly verify this. 
Unlike the case for the unsaturated polyyne radicals, the model results, in agreement 
with observation, show that the the abundance of CH3CCH has the opposite trend 
of lower abundances in lukewarm sources and distinctively larger abundances in hot 
corino sources. The observed abundances of CH3OH show a similar pattern of smaller 
values toward a lukewarm source and larger values toward hotter sources, which the 
model results also accurately replicate. The trends observed in these four species 
provide the framework for a molecular distinction between the lukewarm corino and 
hot corino classes of protostellar envelopes. A recent observation of L1527 indicates 
depletions toward the inner most region, sugg esting that such differentiation occurs 
within an i ndividual source (jSakai et al.ll2010l ). This observation may be consistent 



with that of lOberg et all fl2010h where complex oxygen-containing molecules appear to 
be associated with a cold portion of the low-mass protostar Bl -b. 

The high abundance of C3 observed in absorption in dense regions with T ~ 50 K 
towards W31C and W49N is explained by warm-up models using both our previous 
and updated gas-grain network. Further, this agreement is reached without requir- 
i ng the photo-destructi on of PAH molecules, another suggested source of this species 
(jMookerjea et al.ll2010i ). The new models presented here provide evidence of only a 
slight diminution of abundance as a function of increasing chain length for more com- 
plex bare carbon chains. The formation of C3 from acetylene during the warm-up 
implies the presence of lukewarm corino-like regions within the complex sources. 



Each of these co nclusions cou l d be fu rther tested and improved with more observational 
data. The survey of ISakai et al.l (l2009af ) has indicated that some additional sources such 
as Bl and L1448 N may also cla s sify a s lukewarm corinos on the basis of enhanced C4H 
abundance, while ICordiner et al.l (120111 ) indicated that the enhanced abundance of several 
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carbon chain species toward L1251A (IRS3) could possibly be attributed to a similar effect. 
We encourage follow-up observations of both lukewarm and hot corinos to investigate other 
carbon-chains including CH 3 CCH, as well as CH 3 OH. The depleted abundances of C2H and 
C4H toward hot corinos may make their detection more difficult, but the CH 3 C n H chains 
(n=2, 4, 6, . . .) should be present in enhanced abundance toward these sources. The model 
results for the HIFI sources would benefit from additional detections of longer bare carbon 
chains or an expanded search for other carbon chains in the inventory of lukewarm or hot 
corinos. 

The authors gratefully acknowledge contributions useful discussions of observational 
data with Nami Sakai, Maryvonne Gerin, and Bhaswati Mookerjea. G.H. thanks Rose 
Finn for computational time. E.H. would like to thank the National Science Foundation 
for supporting his research program in astrochemistry and NASA for support to study the 
formation of pre-planetary matter and the analysis of Herschel data. 
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Fig. 1. — Temporal evolution of the C4H abundance for a variety of warm-up scenarios: 
(a) constant T = 10 K; (b) T max = 30 K; (c) T max = 100 K; (d) T max = 200 K. The 
solid lines represent new results including non-zero barrier reactions, while the dashed lines 
represent previous results from lHassel et al.l (120081 ) for the same warm-up conditions. The 
temperature profile for the warm-up is also plotted with respect to time as labeled. The 
observed abundances in L1527 and B228 are indicated with horizontal lines. 
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Fig. 2. — Temporal evolution of the fractional abundances of four species: a.) C2H; b.) 
C 4 H; c.) CH 3 CCH; d.) CH3OH during the warm-up period to T max = 200 K. The solid lines 
represent new results including non-zero barrier reactions, while the dashed lines represent 
previous results from lHassel et al.l ( 120081 ) . The temperature profile for the warm-up is also 
plotted with respect to time as labeled. The observed abundances in the lukewarm corinos 
L1527 and B228 (as I.e.), and the hot corinos IRAS 16293-2422 and for some species NGC 133 
IRAS 4B (as h.c.) are indicated with horizontal lines and grayscale where an observational 
range exists. 




Fig. 3. — Temporal evolution of the abundances of the cyclic and linear forms of C3H for a 
variety of warm-up scenarios: T max = (a) 30 K; (b) 100 K; and (c) 200 K. The solid black lines 
represent the cyclic form, while the dashed black lines distinguish the linear HC3 form. Gray 
long-dashed lines represent the previous results, in which no distinction was made between 
the iso mers. The horizon tal line indicates the observed abundance of the linear isomer in 
L1527 JSakai et al.ll2009bh . 
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Fig. 4. — Temporal evolution of the abundances of the cyclic and carbene forms of C3H2 
for the same warm-up scenarios as in Fig. [3J The solid black lines represent the cyclic form, 
while the dashed black lines distinguish the carbene H 2 C3 form. Gray long-dashed lines again 
represent the previous results, in which no distinction was made between the isomers. The 
horizon tal lines indicate the obse rved abundances of the cyclic and carbene forms toward 
L1527 JSakai et alil2008al . l2009bL as labeled. 
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Fig. 5. — Temporal evolution of the abundances of C5N and HC5N for the same warm- 
up scenarios as Fig. [3j The black solid lines represent the current C 5 N results, the black 
dashed lines represent the current HC 5 N results, and the gray long-dashed and dotted lines 
represent the corresponding HHG results. The horizontal gray band represents the range 
of observations for HC5N toward lukewarm corinos, while the dashed line represents the 
J0rgensen-Bisschop estimate of the upper limit of HC5N toward hot cores. 
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Fig. 6. — Temporal evolution of the fractional abundances of four bare carbon chain species: 
a.) C 3 ; b.) C 4 ; c.) C 5 ; d.) C 6 during the warm-up period to T max = 50 K. The solid lines 
represent new gas-grain results including non-zero barrier reactions, while the dashed lines 
represent models including only gas-phase chemistry. The temperature profile for the warm- 
up is also plotted with respect to time as labeled. The observed abundance of C3 toward 
the massive star forming regions W 31C and W49N is indicated with a horizontal line in the 
first panel (IMookerjea et al.ll2010l ). 
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Table 1. Rate Coefficients for Selected New Classes of Reactions 





Reaction 






n 


a (cm 3 s x ) 


P 


E A (K) 




C n H + H 2 - 


+ C n H 2 - 


h H 


2-9 


1.14(-11) 





950 


2.) 


C„ + H 2 


* CnH 4 


H 


2-9 


1.60(-10) 





1420 


3.) 


C n N + H 2 - 


* HC n N - 


f H 


3,5,7,9 


4.04(-13) 


2.87 


820 
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Table 2. Molecular Fractional Abundances in Lukewarm Corinos 





Observed 




Model 








Species 


L1527 


B228 


T = 10 K 


T max = 30 K 


T max = 100 K 


T max = 200 K 




C 2 H 


1.2(-08)*, a 


5.3(-09) a 


3.0(-10) 


8.9(-09) 


9.9(-09) 


1.0(-08) 


I-C3H 


5.5(-ll) b 




2.4(-10) 


2.3(-10) 


2.5(-10) 


2.6(-10) 


C 4 H 


3.2(-09) a 


2.3(-09) a 


l.O(-lO) 


2.4(-09) 


1.8(-09) 


1.6(-09) 


C 5 H 


1.6(-ll) c 




1.2(-11) 


4-7(-io)i 


3.5(-10) 


3.4(-10) 


C 6 H 


1.0(-ll) c 




1.3(-11) 


2.3(-10) 


1.4(-10) 


1.2(-10) 


C-C3H2 


2.2(-10) b 




6.8(-10) 


4.6 (-09) 


3.8(-09) 


3. 8 (-09) 


H2C3 


1.0(-ll) c 


5.5(-12) a 


4.1(-11) 


2. 5 (-10) 


2.4(-10) 


2. 5 (-10) 


H2C4 


2.7(-ll) a 


6.2(-12) a 


6.5(-12) 


1.6(-10)* 


9.8(-ll) 


8.8(-ll) 


H 2 C 6 


2.5(-12) c 




5.2(-13) 


3.4(-ll) 


1.6(-11) 


1.4(-11) 


CH3CCH 


1.0(-09) a 


1.7(-09) a 


1.2(-10) 


5.3(-09) 


4.3(-09) 


4.3(-09) 


CCS 


8.5(-ll) c 




8.3(-ll) 


l.l(-lO) 


8.1(-11) 


6.8(-ll) 


HC 5 N 


l.l(-10) b 


3.3(-ll) a 


5.0(-13) 


1.7(-10) 


2.4(-10) 


2.4(-10) 


HC7N 


2.7(-ll) c 




4.2(-14) 


3.4(-ll) 


3.8(-ll) 


3.4(-ll) 


HC 9 N 


2.5(-12) c 




2.4(-15) 


6.7(-12) 


5.5(-12) 


4.5(-12) 


HC0 2 + 


1.0(-12) d 


4.5(-12) a 


1.6(-12) 


8.3(-13) 


7.1(-13) 


6.9(-13) 


NH 3 


8.3(-09) e 




2.8(-09) 


4.4(-09) 


6.8(-09) 


1.0(-08) 


CO 


3.9(-05) f 




6.3(-06) 


3.6(-05) 


4.3(-05) 


4.5(-05) 


CN 


8.0(-ll) f 




l.l(-lO) 


l.l(-09) 


2. 7 (-09) 


3.1(-09) 


CS 


3.3(-10) f 




6.1(-10) 


3.3(-09) 


2.4(-09) 


2.0(-09) 


HCN 


1.2(-09) f 




3.2(-10) 


1.2(-09) 


2.8(-09) 


3.8(-09) 


HNC 


3.2(-10) f 




2.8(-10) 


1.1 (-09) 


2.7(-09) 


3. 6 (-09) 


HC3N 


2.0-4.5(-10) b 




9.5(-12) 


5.0(-10) 


9.6(-10) 


1.0(-09) 


HCO+ 


6.0(-10) f 




4.3(-10) 


1.7(-09) 


1.7(-09) 


1.7(-09) 


N 2 H+ 


2.5(-10) f 




5.1(-11) 


l.l(-ll)* 


9.1(-12) 


8.7(-12) 


SO 


1.4(-10) f 




6.6(-10) 


3.2(-10) 


2.5(-10) 


2.1(-10) 


CH3OH 


l.l(-09) b 




3.9(-ll) 


4.2(-10) 


1.8(-10) 


2.1(-10) 


c 3 






6.8(-08) 


6.4(-09) 


6.2(-09) 


5.2(-09) 


*o P t (yr) 






4.2(+04) 


2.7(+05) 


1.8(+05) 


1.6(4-05) 


T (K) 






10 


25 


26 


27 


Kmax 






0.514 


0.599 


0.580 


0.577 


# fits (26) 






19 


19 


20 


19 



*o(-6) = a x lO" 6 

^Results for C4H2 & CeH2 multiplied by 0.02 to approximate abundances of observable carbene isomers 

t Italic type indicates over-production by 1 order of magnitude or more 

tBold type indicates under-production by 1 order of magnitude or more 

feakai et all |2009a| ) 

qSakai et all ll2009bh 

ISakai et all d2008bh 

^Sakai et all l|2008al 'l 

fairota et all l|2009h 

Ijargensen et al.l lj2004h 
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Table 3. Molecular Fractional Abundances in Hot Corinos 





ID AC 1 f}9QQ 9499 


INVjO J-OOO llL/^-O 


NTPr 1 1^^^ TRA^ 
INVjO -LOOO *iLJ 




T = 100 K 


T = 150 K 


T = 200 K 


Species 


^obs 






X 


X 


X 


X 


C 2 H 


5.0(-ll) a 






1.7-2.0(-09p 


2.2-5.6(-10) 


4.0(-12)* 


8.5(-13) 


C 4 H 


3.8(-12) a 




6(-ll) a 


2. 8-3.1 '(-09) 


2.5-4.8(-10) 


4.8(-12) 


1.2(-12) 


CH3CCH 


[2.6±0.3](-07) b (A) 






1.8-24.(-10) 


0.8-1. 2(-06) 


7.4(-07) 


5.0(-07) 


CH3OH 


3(-07) b 


<7(-09) c 




3.7-3.9(-10) 


1.7-8.2(-07) 


2.3(-06) 


1.4(-06) 


CO 


3.3(-05) d 


7.9(-06) d 


1.3(-05) d 


4.3-4.5(-05) 


4.5-4.6(-05) 


4.7(-05) 


4.8(-05) 


CS 


3.0(-09) d 


1.0(-09) d 


1.2(-09) d 


1.2-1.4(-08) 


5.0-5.8(-08) 


5.0(-08) 


5.0(-08) 


SO 


4.4(-09) d 


4.6(-09) d 


3.0(-09) d 


0.1-1. 5(-09) 


4.1-4. 7(- 10) 


1.7(-10) 


2.5(-10) 


HCO+ 


1.4(-09) d 


4.3(-10) d 


6.2(-10) d 


1.3-2. 1(-10) 


0.7-1.9(-10) 


1.6(-12) 


2.0(-12) 


N 2 H+ 


1.4(-10) d 


>1.0(-09) d 


3.2(-09) d 


3.8-4.7(-12) 


3.0-5.0(-12) 


1.3(-13) 


1.7(-13) 


HCN 


l.l(-09) d 


3.6(-10) d 


2.0(-09) d 


0.5-1.0(-06) 


1.3-2.1 (-07) 


3. 9 (-07) 


6.1 (-07) 


HNC 


6.9(-ll) d 


2.8(-ll) d 


1.4(-10) d 


2.5-3.8(-07) 


0.9- 1.5 (-07) 


2. 5 (-07) 


3.1 (-07) 


CN 


8.0(-ll) d 


3.7(-ll) d 


1.4(-10) d 


1.6-2.1 (-07) 


0.3-2. 7(-09) 


1.1(-10) 


6.5(-ll) 


CH 3 CN 


[1.0± 0.4](-09) b 


[1.6± 0.2](-09) c 




4.8-6.3(-09) 


2.1-3.6(-08) 


1.9(-08) 


1.3(-08) 


CH3OCH3 


[2.4± 3.7](-07) b 


<2.8(-08) c 


<1.2-8.0(-08) e 


1.7-4.4(15) 


0.07-6.1(-10) 


1.6(-08) 


1.8(-08) 


CH3CHO 


1.9(-08) b (A) 
3.2(-08) b (E) 






l.l-2.4(-ll) 


0.4-1. 2(-ll) 


8.3(-13) 


7.2(-13) 


H 2 CO 


l(-07) b 


2(-08) f 




l.l-6.0(-08) 


6.8-8.7(-08) 


2.6(-08) 


1.1 (-08) 


HC3N 


1.5(-10) d 


7.2(-ll) d 


l.l(-10) d 


3.4-4.0(-08) 


0.4-1. 2(-07) 


4.9(-08) 


2.8(-08) 


HCOOH 


6.2(-08) b 


[4.6± 7.9](-09) c 




4.2-8. 7(-10) 


0.3-2. l(-08) 


9.6(-08) 


9.2(-08) 


HCOOCH3 


[1.7± 0.7] (-07) (A) b 
[2.3± 0.8] (-07) (E) b 


[3.4± 1.7] (-08) (A) c 
[3.6± 0.7](-08) (E) c 


0.7-3.5(-08) e 
0.5-2. 7(-08) e 


0.8-2. 0(-12) 


0.4-4. 6(-10) 


4.0(-09) 


3.1 (-09) 



HSakai et ail l l2009al) 
HCazaux et"aT1 J2003h 
feottinelli et al.1 ||2004|) 
^Jflrgensen et alj J2004h 
d Sakai et all J2006T) 
iMaret et ail j2004h 

t Italic type indicates over-production by 1 order of magnitude or more 
■fBold type indicates under-production by 1 order of magnitude or more 



